Meson retardation in deuteron photodisintegration above 

7r-threshold t 



O 
=5 



X 



Institut fur Theoretische Physik, Universitat Hannover, D-30167 Hannover, Germany 



M. Schwamb, H. Arenhovel and P. Wilhelm 

Institut fur Kernphysik, Johannes Gutenberg -Universitat, D-55099 Mainz, Germany 

0\ '. Th. Wilbois 

a.- 

(N 

Abstract 

oo 
O 

pled channel approach including iVA- and 7rd-channels with pion retarda- 
0^ ' tion in potentials and exchange currents. A much improved description of 



Photodisintegration of the deuteron above 7r-threshold is studied in a cou- 



total and differential cross sections in the energy region between 7r-threshold 
and 400-450 MeV is achieved. 

PACS numbers: 21.45. +v, 13.40.-f, 25.20.-x 

I. INTRODUCTION 

Despite its long history, photodisintegration of the deuteron is still a very intensely 
studied process. In the past decades, large efforts were made in order to understand 
this simplest photonuclear reaction. At energies below 7r-threshold, where the theory 
is based on realistic iVA-interaction models and corresponding exchange currents, the 
agreement between theory and experiment is quite satisfactory [0-0. On the contrary, 
in the A-resonance region where the A-excitation is dominant, the situation is much 
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less well settled (for a review see 0). None of the models developed so far such as the 
diagrammatic approach of Laget ||, the framework of nuclear isobar configurations in 
the impulse approximation || , the unitary three-body model and the coupled channel 
approach (CC) |7]|| is able to describe in a satisfactory manner the experimental data 
over the whole A-resonance region |3,10]. In the impulse approximation || and also in 



the CC approach of JF|, several parameters were adjusted in order to fit the total cross 
section whereas in the three-body model of [[J and in the coupled channel approach of 
|| all free parameters were fixed in advance by fitting NN- and 7riV-scattering, and tt- 
photoproduction on the nucleon. Consequently, no adjustable parameters remained for 
deuteron photodisintegration. However, it turned out that both approaches considerably 
underestimated the total cross section in the A- region by about 20-30% Another 
failure was the wrong shape of the differential cross section which was flatter than the 
experimental one and which developed a more and more pronounced dip around 90° 
with increasing energy Also the photon asymmetry was not well described, especially at 
photon energies above 300 MeV 0,§. 

One of the crucial points in these calculations is the question of how to fix the jNA- 
coupling Gffif(E/s) in the Ml iVA-current 

j A M N \E A , k) = Gf f^ ] r AN>0 ^ M xfc, (1) 

where E A is the energy available for the internal excitation of the A and k the momen- 
tum of the incoming photon. While in the impulse approximation and also in the 
whole P33 Ml-strength was interpreted as A-excitation strength, the more sophisticated 
approaches |||| have determined G^jv by fitting the M 1+ (3/2)-multipole of pion photo- 
production on the nucleon including nonresonant Born terms, resulting in a smaller G^5v 
and consequently a smaller photodisintegration cross section. The Born terms contribut- 
ing to the (3, 3)-channel are the crossed iV-pole and 7r-pole graphs. When embedded into 
the two-nucleon system, these Born terms become part of the recoil and the 7r-meson cur- 
rent, respectively (Fig. [l]). In static calculations, however, the recoil current is not present 



due to its cancellation against the wave function renormalization current [|TTJ. A similar, 
but less serious problem arises in the treatment of the pion pole diagrams compared to 
the meson current of static MEC. It had already been conjectured in || that this incon- 
sistent treatment of pion exchange may lead to the observed underestimation of the total 
cross section in their coupled channel approach, because by incorporating the Born terms 
effectively into an increased Ml A-excitation strength a satisfactory agreement with the 
data could be achieved. 

In the present paper, we have overcome these shortcomings by including for the first 
time in a coupled channel approach complete retardation in the 7r-exchange contributions 
to potentials and MECs. This retardation is not approximated by keeping only the leading 
order contribution of a p/M-expansion in the 7riViV-propagators as was done for example 
13]. We furthermore consider the intermediate 7r<i-channel in order to ensure all 
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unitarity constraints up to the 27r-threshold. 



II. THEORETICAL FRAMEWORK 

In our framework, the Hilbert space is subdivided into three orthogonal spaces 7i = 
Hjv © 7~C A © 7~Cq HU- T~Ln consists of two nucleons only, in 7i A one of the nucleons is 
replaced by a A(1232)-isobar, and 7{q contains at least one meson besides two baryons. 
TLq is then eliminated by introducing effective operators in TLn®7~La using the projection 
operator technique. 

For photodisintegration we need the outgoing A^A^-scattering wave with total energy 
E which can be expressed as 

|JVJVH) = [l + G%(z)T NN (z) + G*(z)T AN (z)] \NN) , z = E - ie . (2) 

The T-matrix obeys a coupled integral equation of Lippmann-Schwinger type 

T NN (z) = V NN (z) + V NN (z)G% \z)T NN (z) + V NA (z)G£(z)T AN (z) , (3) 
T AN (z) = V AN (z) + V AN (z)G%(z)T NN (z) + V AA (z)G£(z)T AN (z) . (4) 



Here, Gq (z) denotes the free iViV-propagator and Gq(z) the dressed iVA-propagator, 
which depends on the bare A-mass M A and the complex, energy dependent A-self energy 
Xa taken from the dynamical model of Tanabe and Ohta (model A in [|T5|| ) . 

The various effective driving potentials Vxx> with X, X' G {N, A} in Eqs. (|3|) and @ 
consist of two terms (Fig. ||) 

V X x>(z) = V xx ,(z) + VZ d x ,(z) , (5) 

incorporating respectively the usual boson-exchange mechanism and the formation of an 
interacting iViV-pair with deuteron quantum numbers and a pion as a spectator (denoted 
for simplicity by nd). In static applications, we have chosen for V/vat(z) the realistic, energy 
independent Bonn-OBEPR model ||16|| , whereas for the inclusion of full 7r-retardation we 



have used an improved version of the energy dependent Bonn-OBEPT model developed 
by Elster et al. jnj which contains besides retarded operators also self energy diagrams 



calculated within a simple Lee model |L8 



For V/va and Vatv = Vnai we have constructed corresponding static and retarded 7r- 
and p-exchange potentials, taking for the kNN- and 7riVA-form factors the ones of the full 
Bonn potential ]l6j. The pNN- and piVA-regulator masses are fixed by fitting the x D<i 
iVTV-partial wave which is of crucial importance for deuteron photodisintegration because 
of its strong coupling to the dominant 5 S , 2(iVA) partial wave. The diagonal interaction 
Vaa i n consists of the 'forward' and 'backward' going pion diagrams as depicted in 
Fig. |2|. For its 7riVA-vertex, the form factor of [0 is used. 

We now turn to the discussion of Vx X ,, which was studied in but not included in 
M. From the corresponding diagram of Fig. [2] it is obvious that V xx , can be written as 

Vxx'i z ) = Vx'qG t QQ (z)G V Q x> , (6) 

where V x >q is the nonrelativistic NNtc- or iVA7r-vertex, respectively. The free irNN- 
propagator Q^ NN > stat j s taken in the static limit and t 1 QQ{z) denotes the iViV-scattering 
amplitude in Q-space describing an intermediate off-shell iViV-state with the quantum 



numbers of the deuteron and a pion as spectator. The amplitude t^n(z) is obtained from 
the Lippmann-Schwinger equation 

tl%(z)=v* d + v« d Gf(z)F Q d Q (z) , (7) 

where Gq 6, denotes the free 7rc?-propagator. It can be solved easily by assuming a separable 
form for the driving term v %d 



d _V NN {E^)\d,-q){d,-q\V NN {E^) ^ _ _ ^ ^ q ^ _ ^ _ 
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where |d7r;<f ) denotes a free 7rc?-state with relative momentum q. This ansatz for v nd 
satisfies for any realistic iViV-potential Vnn the Schrodinger equation for the Tid-state 

G% d (E nd )v wd \dTr; q > = \dii; q) . (9) 

Since the parameters of the realistic Vnn are fitted to deuteron properties and NN- 
scattering phases without explicit coupling to the NA- and ird- channels, we have to apply 
a box-renormalization 

Vnn(z) Vnn(z) - Vna(2M)G£(2M)V an (2M) - V^ d N (2M) (10) 

in order to ensure an approximate phase equivalence in the presence of these additional 
channels. 

Now we turn to the e.m. part of our model. Above 7r-threshold, the A-excitation 
is the most important photoabsorption mechanism, which is described by the current 
operator in Eq. ([I]) neglecting small E2 contributions. As mentioned above, G^ N (E A ) is 
an energy dependent and complex coupling which is fitted to the experimental M 1+ (|)- 
multipole of pion photoproduction on the nucleon. It contains besides the bare ^NA- 
coupling contributions from nonresonant pion rescattering. The full pion photoproduction 
amplitude t ni (E A ) in the (3, 3)-channel can be written as 

U,(E A ) - t W7 (E A ) - £ a _ m o_ Sa(Ea) , (11) 



where t% L(-Ea) is the nonresonant Born amplitude. While in || an effective 7iVA-coupling 
^ajv(-^a) and the model of [14] for the bare A- mass M^, the A-self energy Sa and the 
A7riV-vertex v* A has been used, we follow here the work of Tanabe and Ohta (model A in 



||15|| ) in which an explicit evaluation of the half-off-shell Born amplitude in the M 1+ (3/2)- 
multipole has been performed. 

Furthermore, concomitant with the construction of the effective interactions, we have 
derived the corresponding effective two-body charge and current density operators which 
fulfil current conservation (for details, see These 7r-retarded MECs contain besides 

the usual vertex-, meson- and contact-MECs the recoil current and charge densities and 
a couple of additional two-body operators which vanish identically in the static limit. 
Whereas the latter ones yield only very small contributions which can be safely neglected, 
the recoil contributions turn out to be quite important (see discussion below). They 
do not appear in static approaches due to their cancellation against the wave function 
renormalization contributions |Tl|, which have their origin in the renormalization of the 
baryonic states when eliminating the mesonic wave function components. This concept 
breaks down beyond the 7r-threshold if full 7r-retardation is considered since the nNN- 
component can be on-shell. Therefore, we do not orthonormalize and no wave function 
renormalization contributions appear. Consequently, the recoil current and charge densi- 
ties have to be included. 

Besides the 7r-MECs and the corresponding Siegert operators, our model includes the 
usual nucleonic one-body current and as most important relativistic contribution (besides 
retardation) the spin-orbit current. Because the p-mass is rather large, retardation in the 
p-MEC is expected to be rather unimportant and therefore not considered in this work. 
Concerning the iVA-7r-MEC, we include retardation in the corresponding recoil part, 
whereas, due to their minor importance, the contact-, meson- and vertex-contributions 
are taken in the static limit. Since the pion production model of Tanabe and Ohta [|i~5|l 
effectively incorporates tu-exchange, we include in addition the leading order pir^- and 
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co>7T7- currents, which are purely transverse 



go 



III. RESULTS 

In the numerical evaluation we have included all multipoles and scattering waves 
up to L = 4 and j = 4. Retardation in the meson-, vertex-, contact- and pn^juj-n^- 
currents is only included up to L = 2. For the higher multipoles the static limit has been 
adopted, because inclusion of retardation effects there does not show significant effects 
on observables. We start the discussion of deuteron photodisintegration with the total 
cross section shown in Fig. |3]. Similar to |8|||, the static calculation (without 7r<i-channel 
and p7T7/c<j7r7-MECs) considerably underestimates the data (dotted curve). However, the 
7iVA-coupling of used in our approach is somewhat larger than the one used by 
Wilhelm et al., so that the discrepancy between our static calculation and experiment 
is not as dramatic as in The reduction of the cross section by the inclusion of 

retardation in the hadronic interaction (dash-dotted curve) is, however, overcompensated 
if retardation in the 7r-MECs is added, of which the recoil contributions turn out to be the 
most important ones. This destructive interference of retardation in potential and recoil 
current corresponds to the cancellation of wave function renormalization and recoil current 
in the static limit. However, here remains a net effect leading to a significant increase. 
The cross section is further enhanced by inclusion of the 7r<i-channel and the p-K^/uj-K^j- 
MECs (by about 8% respectively 3% at 260 MeV) so that the complete calculation (full 
curve) agrees quite well with the experimental data over the whole energy range. 

Concerning the role of the 7r<i-channel in iViV-scattering, its influence on the inelas- 
ticity of the 1 D 2 -channel is similar to ||. However, in contrast to our result of a slight 
enhancement, the 7r<i-channel leads in [|J to a reduction of the cross section in 7c? — > pn, 
the reason for which we are presently investigating. 

In Figs. f| and |^, we show differential cross sections and photon asymmetries for various 
energies. Whereas now the differential cross section is in satisfactory agreement with the 
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data, we slightly underestimate the absolute size of the asymmetry. In contrast to HH we 
are able to reproduce quite well the shape of these two observables also at higher energies, 
in particular the dips in the angular distributions have disappeared. 

In conclusion, we have demonstrated that a CC approach including NA- and ird- 
channels with full pion retardation in potentials and exchange currents is able to remove 
almost quantitatively the still existing discrepancies between theory and experiment in 
deuteron photodisintegration in the A-region. Further improvements of our model should 
include relativistic contributions to internal dynamics and currents whereas boost effects 



are expected to be small [20]. Furthermore, a realistic NN- interaction is desirable which 
includes isobar degrees of freedom in a coupled channel approach from the beginning so 
that no box renormalization procedure is necessary. The success of our model in —>■ NN 
encourages us to study retardation in related e.m. processes on the deuteron in the near 
future, like 7r-photoproduction and electrodisintegration. 
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FIGURES 
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FIG. 1. The Born terms contributing to the Mi + (|)-multipole amplitude of pion photo- 
production (upper part) and their correspondence in the recoil and pion-MECs in deuteron 
photodisintegration (lower part). 





X' 


N 




X' 


X' 


N 












7t 






N 


X 




N 


X 


N 



NNpS^'D,) 



FIG. 2. Graphical representation of the hadronic interaction of Eq. (||) with X, X' € {N, A}. 
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FIG. 3. Total cross section for 7 a! — > pn as a function of photon energy in compari- 
son with experiment ||.1C|. Dashed: static calculation of Wilhelm et al. Q: dotted: static 
OBEPR-calculation in our approach; dash-dot: retardation switched on in the hadronic part 
only, but static MECs; full: complete calculation including 7rci-channel and /37r7/u;7r7-MECs. 
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FIG. 4. Differential cross sections for various energies in comparison with experiment f9|,|lQ] . 
Notation of the curves as in Fig. 0. 
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FIG. 5. Photon asymmetry £ for various energies in comparison with experiment 0,^1,22] 
Notation of the curves as in Fig. 0. 
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